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ABSTRACT 

The Las Campanas Redshift Survey (LCRS) contains 23697 galaxies, with an 
average redshift z = 0.1, distributed over six 1.5° by 80° slices in the North and 
South galactic caps. We have computed the power spectrum P(k) for magnitude- 
limited samples of LCRS galaxies over wavelengths A = 2n/k = 5-400 hr 1 Mpc. 
The LCRS P{k) may be approximated as oc k~ 18±0A for small scales A = 5 — 
30 h~ x Mpc, changing to oc k x±1 for large scales A ~ 200 — 400 h~ x Mpc. The 
overall amplitude corresponds to er 8 = 1.0 ± 0.1 in redshift space. 

Comparisons to the power spectra of other redshift surveys will be pre- 
sented; the LCRS results agree best with those from the combined Center for 
Astrophysics (CfA2) and Southern Sky redshift surveys (SSRS2). For A ^ 
100 h~ x Mpc, the LCRS results are consistent with those of other surveys, given 
the large errors among all the surveys on these scales. For A <^ 100 h^ 1 Mpc, the 
LCRS P(k) is well determined and similar in shape to the P{k) of other surveys, 
but with an amplitude differing from some of the other samples, possibly because 
of inherent clustering differences among different types of galaxies. In particu- 
lar, power spectrum measurements for volume-limited LCRS samples show that 
galaxies brighter than about M* — 1 appear about 50% more strongly clustered 
than those fainter. Also, a sample of LCRS emission galaxies shows 30% weaker 
clustering than the full LCRS sample. 

Comparisons to N-body models show that the LCRS power spectrum lies 
intermediate between that of a standard flat f2 /i = 0.5 cold dark matter (CDM) 
model and an open QqIi = 0.2 model, both normalized to ag — 1 for galaxies. 
On large scales A ^ 40 h~ x Mpc, we have fit the LCRS results to various linear 
CDM models, and find that a number of them could meet the constraints set 
by the LCRS power spectrum, the Hubble constant range 0.5 <^ h ^ 0.8, the 
abundance of galaxy clusters, and the reasonable assumption that LCRS galaxies 
are roughly unbiased tracers of the mass, relative to the normalization provided 
by the 4-year COBE DMR data. The possibilites include open CDM or flat non- 
zero cosmological-constant CDM models with Q « 0.4—0.6 and shape parameter 
T Q h ~ 0.2 — 0.3, as well as flat Qo = 1 models with massive neutrino density 
Vt v ~ 0.2 - 0.3 or a spectral tilt n ^ 0.7 - 0.8. 

Subject headings: cosmology: observations — galaxies: clustering — large-scale 
structure of universe 
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1. Introduction 

The power spectrum of density fluctuations is an important fundamental quantity of 
interest for the problem of structure formation in cosmology. For example, for Gaussian 
density fluctuations, as may arise out of inflationary scenarios for the origin of the universe 
( |Bardeen, Steinhardt, fc Turner 1983|) , the power spectrum provides a complete descrip- 
tion of the initial fluctuations. The subsequent development and evolution of structure via 
gravitational instability, in the context of a dark matter dominated universe and models of 
biased galaxy formation QBardeen et al. 1986|) , have received much theoretical attention in 



the past decade, both analytically and through the use of large numerical N-body simula- 
tions (e.g., Davis et al. 1985 ; Park 1991 ; Zurek et al. 1994 ). On the observational front, 



the detection of microwave background anisotropies by the COBE DMR experiment ( |Smoot 



et al. 19§2] ) provided important constraints on the primordial power spectrum on horizon 



scales. The clustering of galaxies on smaller scales and the interpretation of the observations 
to reveal the matter content of the universe, to elucidate the process of structure formation, 
and to connect to the primordial fluctuations seen by COBE, remain important problems in 
observational cosmology (see e.g. [Lyth fc Liddle 1994| for a review). 

Much recent attention has been paid toward determination of the power spectrum of 
galaxy clustering in a variety of galaxy surveys, including the CfA ( |Vogeley et al. 19"§2"| ; [Park] 



|et al. 1994[ ), SSRS ( [Park et al. 195% |da Costa et al. l994b| ), IRAS 1.2 Jy ( [Fisher et al. Iggg ), 
IRAS QDOT ([Feldman et al. 1994ft , and APM ([Baugh fc Efstathiou 1993; , 1994) surveys. In 
brief, the power spectra of these surveys have appeared inconsistent with predictions of the 
"standard" biased cold dark matter (CDM) model of structure formation with Q h = 0.5 
( piumenthal et al. 1984|) , while an unbiased QqIi ~ 0.2 model with more large scale power 



agrees better with the observations (e.g. |da Costa et al. 1994b| ). (We express the Hubble 



constant as Hq = 100 h km s _1 Mpc -1 , and will use h = 1 unless otherwise indicated.) In 
this paper we present the power spectrum for galaxy samples drawn from the Las Campanas 
Redshift Survey (LCRS), an optically-selected survey of 23697 galaxies with an average 
redshift z = 0.1. The large sample size and extent of our survey allow us to examine 
the power spectrum up to wavelengths of ~ 400 h~ x Mpc, and to provide measurements 
independent of previous results for the purpose of comparing against cosmological models. 
In particular, measurements of the power spectrum on the largest scales A ^ 100 h^ 1 Mpc are 
especially interesting, as we expect the power spectrum to peak there and begin its turnover 
toward the primordial spectrum constrained by COBE and other microwave background 
observations. The precise amplitude and shape of the power spectrum on large scales will 
provide important clues in discriminating among cosmological models. 

A detailed description of the Las Campanas survey is given in Shectman et al. (1996), 
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and additional particulars may be found in Shectman et al. (1992, 1995), Tucker (1994), Lin 
et al. (1996) and Oemler et al. (1993). Here we briefly describe the main survey parameters. 
The survey geometry is that of six 1.5° x 80° "slices" (declination by right ascension), three 
each in the North and South galactic caps. Figure |I] shows the LCRS galaxy distribution and 
clearly illustrates the striking pattern of clusters, filaments, walls and voids that is present. 
The first 20% of the data was obtained using a 50-object fiber-optic spectrograph, and the 
remaining 80% of the data was taken with a 112-object system. The nominal isophotal 
magnitude limits for the 50-fiber data were 16.0 < m < 17.3 ("hybrid" Kron-Cousins R 
magnitudes), and an additional cut was applied which excluded the lowest 20% of galaxies 
by central surface brightness. For the 112-fiber data, the nominal magnitude limits were 
15.0 < m < 17.7, with exclusion of just the lowest 4-9% of galaxies by surface brightness. 
The survey photometric limits were chosen so that there would be typically more targets 
per field than available fibers, and we selected targets at random among those that met the 
selection criteria. The survey slices were built up by observing 1.5° x 1.5° fields, one at a time, 
with a maximum of 50 or 112 galaxies observed per field. Because we generally do not re- 
observe any of our fields, we must keep track of the variable field-to-field sampling fractions / 
in our subsequent statistical analyses. The average sampling fraction is 70% for the 112-fiber 
data and 58% for the 50-fiber data. Also, mechanical constraints prevent two object fibers 
in a single spectroscopic field from approaching closer than 55", introducing an additional 
geometric selection effect. We will find below that the various sampling, photometric, and 
geometric selection effects in our survey do not significantly affect the power spectrum results. 

In § ^| we detail our power spectrum estimation techniques and verify them on N- 
body simulations. In § ^ we present the power spectra of magnitude-limited samples of Las 
Campanas galaxies, and compare our results to the power spectra derived from other redshift 
surveys. In § [| we compute the power spectrum for volume-limited samples of Las Campanas 
galaxies and test for luminosity bias in the survey. In § [5] we compare our power spectrum 
results against those from N-body simulations. We will then focus on the large scale linear 
power spectrum, relate our results to the COBE DMR constraints, and compare against 
the predictions of several classes of CDM models. We summarize our results in § ^. Note 
that a complementary analysis of the 2- dimensional LCRS power spectrum has already been 
carried out QLandy et al. 1996| ; more on this below), and the derivation of the closely- related 
two-point correlation function of LCRS galaxies is described in Tucker (1994) and Tucker et 
al. (1996). 
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2. Estimating the Power Spectrum 

The power spectrum estimation technique used here has been described by various 
authors, in particular see Fisher et al. (1993), Feldman et al. (1994), and Park et al. (1994), 
and we detail the method below. The most important difference is that the LCRS consists of 
six essentially two-dimensional "slices" , so that we need to account for "convolution" effects 
caused by the survey geometry in order to calculate the power spectrum properly. These 
convolution effects are also evaluated below. 



2.1. Methods 

Given a galaxy density field p(r) with density contrast 5(r) = [p(r) — p]/p, where p is 
mean galaxy density, we have the two-point correlation function 

Z(r) = (6(*'W + r))r> (1) 

and the power spectrum 

P(k) = J d 3 r £(r)e ik - r . (2) 

With this definition P(k) is a power spectral density with units of volume. Assuming isotropy 
we have f (r) = £(r = |r|) and P(k) = P(k = |k|). 

Now, suppose we have a survey of volume V, with N observed galaxies at positions rj. 
The survey's selection function s(r) is defined to be the fraction of galaxies at position r 
expected to be observable by the survey, given the survey's photometric and other selection 
criteria. Also, let w(r) = l/s(r). We can then form an estimator <5(r) for the true density 
contrast 5(r) by 

' (r) ^E»F| rf,(r - r ' ) - 1 - <3) 
where 6^ is the Dirac delta function and where we have estimated the mean density p by 
E?=i w ( r i)/ V - Note that 

we will use a circumflex " to denote quantities estimated from the 
observations. We next take the Fourier transform 

8(k) = l/rf 3 r5(rK kr (4) 

1 N 

Y^wir^-Wik) , (5) 



where W(k) is the survey window function 

W(k) = ^J v d 3 r e ik ' r . (6) 
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We can convert the sum over galaxies in equation (|5|) to a sum over infinitesimal cells 
distributed throughout the survey volume, as in Peebles (1980, §§ 36 and 41), where each 
cell i holds rij = 1 or galaxies. Then the expectation value of |<5(k)| 2 may be written 

(I WD = ppEEW^W^^e*'^-^ (7) 

+ l^(k)l 2 

r i 



Using equations (0) and (|), and the relations 

(mnj) = 
we can rewrite equation (|7[) as 



p 2 [l + £(|r, - T^MrMr^dVidVj % ± j 
ps{ri)dVi i = j 



<iwr> = ^ + (s) , (9) 

which is the sum of the true power spectrum convolved with the survey window function, 

P(k) EE / d 3 k' P(k')\W (k-k')l 2 , (10) 



(2tt) 3 „ 

and the expectation value of the shot noise arising from sampling a finite sample of galaxies, 

(S) = ^iJ d 3 r Mr) . (11) 

Note that we are using a tilde ~ to denote the convolved power spectrum. For a given set of 
galaxies the shot noise S is given exactly by the i = j terms in the first sum of equation (0), 

^IeWS^' (12) 

and our first estimate of the observed power spectrum, convolved with the survey window 
function, is thus 

dQk 



P(k) = Vj^[\5(k)\'-S], (13) 

where we average over different wavevectors k at fixed magnitude |k| = k. However, because 
we do not know the true mean density beforehand, but must estimate it from the survey itself, 
we will underestimate the power spectrum on scales comparable to the survey size by an 
amount P(0)|TU(k)| 2 ( [Peacock fc Nicholson 19911 ). It turns out that for our survey geometry, 
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P(k) ~ constant « P(0) on large scales (see the next subsection). Given an initial estimate 
for P(0), we can correct for the underestimate by adding back the P(0)|W / (k)| 2 term, so 
that our final estimate for the observed convolved power spectrum is 

Hk) = / ^ (V [|5(k)| 2 - S] + P(0)|jy(k)| 2 ) . (14) 

We limit our power spectrum computations to scales for which this correction is <^ 10%, 

which corresponds to A ^ 400 h^ 1 Mpc for the LCRS. Now, the expectation value of P(k) 
is, using equation fllCf), 

(Hk)) = (15) 



where 



dk' k' 2 P(k')K(k,k') , (16) 

^,fcO--^/^/^|^(k-kO| 2 . (17) 

The effect of the convolution integral ( |16"D for our survey will be illustrated in the next 
subsection. Because convolution effects are large for our survey, it will be convenient to 
deconvolve the observed P(k) to recover P(k). We will use the iterative method due to 
Lucy (1974), which has been applied in a similar context by Baugh and Efstathiou (1993) 
to recover the spatial power spectrum from the angular correlation function. Given the 
observed convolved power P at a set of wavenumbers kj, we can find the true power P at the 
wavenumbers ki by starting with initial guesses P°(ki), and then computing new estimates 
P n (ki) by iteration 

where 

P n (kj) = E k 2 l P n {k l )K{k v k)Ak . (19) 
i 

We also apply a smoothing procedure 

P n (h) -> O.SP^) + 0.25 [P n (fcj_!) + P n (k i+1 )] (20) 

in order to reduce noise as the P n (ki) converge to a final solution. On the largest scales 
k ^ 0.03 h Mpc -1 , the smoothing also reduces the sensitivity of the solution to the choice of 
initial guess P (fc») (we adopt P°(fcj) = constant). Because the deconvolution procedure is 
not too well constrained on large scales, and because the smoothing procedure is somewhat 
arbitrary, we will use the deconvolved power spectra only for qualitative comparisons. The 
computation of the convolved power spectrum from a given true power spectrum is more 
straightforward, and our quantitative analyses will focus on the convolved power spectrum. 
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2.2. Test of the Methods on N-Body Simulations 

We check our power spectrum estimation methods on an N-body simulation kindly 
provided by Changbom Park. The simulation is an open cold dark matter model (denoted 
ODM), with Q = 0.4 and h = 0.5, whose power spectrum has been found to be a good 
match to that of the Center for Astrophysics (CfA) Redshift Survey (|Park et al. 1994f ). The 



model is unbiased, that is, b = 1, where b is the ratio of galaxy to mass density fluctuations. 
The model is normalized so that erg = 1, where os is the rms galaxy density fluctuation in a 
sphere of radius 8 h^ 1 Mpc. The ODM simulation contains 240 3 particles and was computed 
using a particle-mesh code ( |Hockney fc Eastwood 1981| ) on a 480 3 mesh with corresponding 



physical comoving volume (576 h 1 Mpc) 3 . 

We draw from the ODM simulation a set of "mock" redshift surveys using the geometric 
and photometric selection criteria of our actual survey (see § |3.1|) . The ODM "galaxies" are 



assigned absolute magnitudes using the luminosity function derived for the LCRS ([Lin et al. 



1996|) . The velocities of the ODM galaxies are included so we measure the power spectrum in 



redshift space. For each mock survey we rederive the selection function and then compute the 



power spectrum using the methods detailed in § [2.1[ Now, even though the ODM simulation 
box is quite large, the full LCRS still only fits along a body diagonal of the ODM box. We 
can however fit either the Northern or Southern halves of the survey into the box fairly 
readily, and we have drawn a series of Northern and Southern LCRS mock surveys. To make 
a reasonably large number of whole-LCKS mock surveys, we chose to combine a Northern 
and a Southern mock survey even though these two mock surveys are not originally oriented 
relative to each other in the simulation box as the real Northern and Southern surveys are. 
Figure |2| shows the average power spectrum measured from 30 whole-LCRS mock surveys. 
We show both the directly measured, or convolved, power, as well as the deconvolved power 
of the mock surveys. The errors plotted are the standard deviations of the mean (sdom) 
of the 30 surveys. Also plotted in the figure are the true unconvolved power spectrum of 
the ODM simulation and the ODM power spectrum after convolution with our whole-survey 
window function (using equations |]TB[ and |T7[). There is good agreement between the true 
ODM power spectrum and the average power spectrum of the mock surveys over the range 
of scales A = 5 - 400 h' 1 Mpc that we examine, showing that we can correctly account for 
the geometry and selection criteria of the LCRS in our measurement and deconvolution of 
the mock-survey power spectra. 

Note from Figure |^ that the main effect of the convolution by the window function of 
our survey is that the convolved power spectrum lies below the unconvolved power spectrum 
over most of the scales probed. Also, note that on the largest scales, wavelengths A > 
100 h~ l Mpc, the convolved power spectrum is flat whereas the true power spectrum may 
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peak and begin to turn over. (This is seen more clearly for the power spectrum of Figure |6|.) 
In Figure ^ we illustrate the convolution effect in more detail by plotting the integrand 
k' 2 P(k')K(k, k') of the convolution integral ( |TE| ) for the ODM model at several values of k. 
This integrand shows the contribution to the convolved power P at k from the true power 
P at different k! . For a 3-dimensional survey geometry, such as a sphere or cube, at scales 
A = 2ii/k small compared to the survey size, we would expect this integrand to be sharply 
peaked at k' = k\ that is, convolution effects would not be important and the measured 
convolved power would be very close to the true power. However, the Las Campanas survey 
has essentially a two-dimensional geometry: it is large in two dimensions but thin along the 
third (cf. Figure [I]). Convolution effects can be noticeably important at scales smaller than 
the survey size, and "aliasing" occurs: power at k' ^ k makes a significant contribution to 
the convolved power P at k, leading to the effects noticed from Figure fj. As seen in Figure [5], 
for k ^ 0.1 h Mpc -1 , the integrand is broad and somewhat complicated in shape, including 
contributions from a large range of k! . The situation improves for smaller scales and larger 
k, as the integrand becomes more and more sharply peaked at k' = k. Nevertheless, though 
these convolution effects are complicated in detail, they can be computed and accounted for 
in the power spectrum analysis, as we saw earlier. 



3. Power Spectrum of LCRS Galaxies 

3.1. Power Spectrum of LCRS Samples 

We compute the power spectrum in redshift space, but first convert our heliocentric 
redshifts to comoving distances r, assuming Qq — 1. Using a different value for Qq does not 
make much difference; for example, r changes by 2% at z = 0.1 if we instead use Qq = 0.2. 
We first consider the power spectrum for magnitude-limited samples of the LCRS. We thus 
need to weight by the inverse of the selection function, but we note that this weighting 
becomes large when the selection function becomes small at large distances, where only 
the intrinsically brightest and rarest galaxies are observed. Also, because of the surface 
brightness cuts imposed on the sample, we do not probe the galaxy luminosity function 
well at the faint end ( |Lin et al. 1996|) , but these intrinsically faint galaxies could make a 
large contribution in the weighting. To reduce the sensitivity of our power spectrum results 
to these potential problems, we make additional restrictions on the redshifts and absolute 
magnitudes of the galaxies used in our clustering analyses: 10000 km s _1 < cz < 45000 
km s' 1 (or 98 < r < 405 h" 1 Mpc) and M 1 = -22.5 < M - 5 log h < -18.5 = M 2 . These 
limits keep the selection function s(cz) ^ 0.05, and result in a total sample of 19305 galaxies. 
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The weight function w used in equation (|5|) for galaxy i of the survey is given by 

»<«> = ^ • (21) 

where is a weighting factor which accounts for variable field-to-field sampling fractions 
and incompletenesses as a function of apparent magnitude and central surface brightness 
(details in [Lin et al. 1996Q , and s(r$) is the selection function appropriate to the survey field 
containing galaxy i, 

r min[M2(ri),M 2 ] I r M 2 

sin) = <f>(M)dM / / <f>(M)dM . (22) 



max[Mi(ri),Mi] / J Mi 

Here Mi(rj) denotes the minimum and M 2 (rj) the maximum absolute magnitude observable 
at the distance r, given the apparent magnitude limits of the field of galaxy i, and 4>(M) is 
the LCRS luminosity function. For 4>(M) we use a Schechter luminosity function ( [Schechterj 
1976 ) with parameters M* = —20.3 + 5 log h and a = —0.7; details concerning the derivation 



of the luminosity function are given in Lin et al. (1996). 

Figure |](a) shows the observed convolved P(k) for the full LCRS sample defined above, 
computed using the methods detailed in § |. We have also split the full sample into four 
subsamples by dividing each of the Northern and Southern data subsets in half by RA; see 
Table [3] for details. The individual P(k) for the four subsamples are also plotted in Figure [|. 
We use the standard deviations of the mean of P(k) for the subsamples, or the la errors from 
the whole-LCRS ODM mock surveys of § |2.2| , whichever is larger, to estimate the la errors 
for the full sample. We also note that though we average P(k) over different directions of |k| 
at fixed |k| = k in equation (0), results from the mock surveys show that the errors depend 
only weakly on direction. The convolution effects are quite similar for all five samples shown 
in Figure £|(a), so that the convolved power spectra may be directly compared. There is good 
agreement among the four subsamples for wavelengths A ^ 100 h~ l Mpc. On larger scales, 
aside from sampling fluctuations, uncertainties in the survey's selection function and mean 
density will also contribute to the differences among the samples. We have also divided 
the data into 50-fiber and 112-fiber subsets (see Table [I]), and computed the convolved 
power spectrum for each subsample, with the results plotted in Figure |](b). Again, the 
convolution effects for these two subsets are similar so we directly compare them. There 
are no conspicuous systematic differences between the power spectra computed from 50- vs. 
112-fiber data, giving us confidence that the power spectrum is not sensitive to the different 
sampling and photometric selection criteria involved. The two subsamples agree well for 
A £ 150 h- 1 Mpc. 

We can also check the consequences on the power spectrum arising from two geometric 
selection effects of the survey. First, recall from § |l]that an instrumental constraint prevents 
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two galaxy fibers from approaching closer than 55" on the sky. In the ODM mock surveys we 
have actually applied this restriction and our previous results indicate apparently little effect 
on the power spectrum. More explicitly, we can assign fake velocities to those unobserved 
objects which lie within 55" of an observed survey galaxy, and then recompute the power 
spectrum. Some fraction of these pairs will be physically associated, and we assign velocities 
to the unobserved objects according to the distribution of pairwise velocity differences for 
observed survey galaxy pairs with angular separations between 1' and 2'. For this purpose we 
use a very generous upper bound of 5000 km s^ 1 for the velocity difference, corresponding to 
60% of the observed velocity difference distribution. The remaining 40% of pairs are assumed 
to be chance superpositions, for which we assign velocities to the unobserved objects accord- 
ing to the overall redshift histogram of the particular survey field containing that pair. Note 
that this procedure adds an average of 1351 galaxies to the full LCRS sample used for power 
spectrum computation. Figure ||(a) plots the full range of the power spectrum computed 
from 10 different sets of such velocity assignments, and confirms that the full-sample LCRS 
power spectrum is indeed very little affected over the scales probed (< 10% difference). Sec- 
ond, as described in Shectman et al. (1995), the spectroscopic success rate is lower (though 
still > 90%) at the corners of the spectroscopic fields. We can approximately compensate for 
this effect by calculating a "neighborhood" sampling fraction within an area of radius 1000" 



centered on each galaxy, and then applying this additional geometric weighting in fl2~T|) . The 
choice of regions of 1000" radius to define the weights came from trial and error efforts to best 
reproduce the angular pair distribution of the survey photometric catalog, by appropriately 
weighting the pair distribution of the spectroscopic sample. Figure |5|(b) shows that applying 
this correction also has little effect on the LCRS power spectrum. For simplicity, since both 
the above geometric effects are small, we will neglect them in the remainder of this paper. 

We next use the Lucy deconvolution procedures outlined in § ^ to obtain estimates of 
the deconvolved P(k) for the full sample and the four subsamples, with the results for the 
full sample displayed in Figure Here we again estimate errors using the sdom of the four 
subsamples, or the la errors from the deconvolved power spectra of the whole-LCRS ODM 
mock surveys, whichever is larger. Besides using Lucy deconvolution, we have also convolved 
a convenient model for the power spectrum and fit to the observed LCRS full sample P(k): 

P(k) ~ 21X2 {k/ko)n ° (23) 

The model was motivated by a similar fitting formula applied to the APM data by Peacock 
(1991) (see also [Feldman et al. 1994|) . The model has four free parameters (ko,n ,k c ,n c ), 



and essentially it behaves as a power law with slope n — 3 at small scales (large k), changing 
to another power law with slope n — 3 + n c at large scales (small k) (we take no,n c > 0). 
The transition occurs near wavenumber k c , and there remains an overall normalization that 
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can be more meaningfully expressed in terms of erg. We find 

k c = (0.06 ±0.01) h Mpc" 1 

P(k) oc k n °- 3 = A;- L8±ai , 27i/k = 5-30 hT l Mpc 

P(k) oc k no ~ 3+nc = k x±1 , 2-n/k « 200 - 400 h" 1 Mpc (24) 

k = (0.17 ±0.01) h Mpc^ 1 

cx 8 = 1.0 ± 0.1 (redshift space) 

We plot the fit in Figure ||. Note that for A £ 100 hr x Mpc, the power spectrum is well- 
determined. There is an apparent peak at about 2ir/k c = 100 h~ l Mpc, but on larger scales 
the form of P(k) is not too well determined, with our fit power law slope consistent with a 
flat to declining P{k) with decreasing k. We want to emphasize that for A ^ 100 h^ 1 Mpc, 
the errors are unavoidably large and we do not claim detection of a turnover in the power 
spectrum, only that our results suggest that P(k) is at least flattening out on these largest 
scales. Improved power spectrum estimation techniques, as recently described by Tegmark 
(1995) or by Vogeley & Szalay (1996), may prove to be helpful in reducing the errors on the 
large-scale power spectrum, though application of these methods is beyond the scope of the 
present paper. 

As we mentioned in § [lj, a complementary analysis of the 2D LCRS power spectrum 
may be found in Landy et al. (1996). The wavevectors in the 2D treatment are aligned along 
the "planes" of each of the six slices, rather than averaged over the whole sphere as they 
are in this paper. The 2D analysis thus emphasizes those spatial directions for which the 
survey geometry is best suited on large scales. This renders the 2D analysis less sensitive to 
the presence of aliased power and therefore more sensitive to the presence of features in the 
power spectrum on large scales. In particular, a striking peak at A ~ 100 h~ l Mpc is seen in 
the 2D power spectrum, and it is a signature of the clear pattern of over- and under-dense 
regions (e.g., walls and voids) in the LCRS galaxy distribution (Figure [I]). A different "core- 
sampling" analysis applied to the LCRS by Doroshkevich et al. (1996) also indicates that 
100 h^ 1 Mpc is approximately the typical separation between large sheet-like structures in 
the galaxy distribution. These results are also reminiscent of a similar peak at 128 h~ l Mpc 
in the ID power spectrum of deep pencil beam surveys analyzed by Broadhurst et al. (1990). 
A peak at A ~ 100 h~ l Mpc is also suggested in the 3D power spectrum fit above, though it 
does not stand out in Figure |6], both because of the logarithmic y-axis scaling and because 
of the suppression due to aliasing as mentioned above. The focus of the 2D paper is on 
the 100 h- 1 Mpc scale power, and the reader is referred to Landy et al. (1996) for the full 
details. On the other hand, the current analysis has the broader goal of relating the full 
LCRS power spectrum, in the more usual 3D context, to results from previous surveys, to 
COBE constraints, and to predictions of cosmological models. 
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3.2. Comparison to Other Surveys 

We proceed now to compare the LCRS power spectrum to that determined from several 
other large galaxy redshift surveys. First, Figure |7| shows the comparison to the power spec- 
tra of three samples drawn from the SSRS2 and CfA2 redshift surveys: two volume-limited 
samples SSRS2+CfA2 101 (r < 101.05 h~ x Mpc, M B < -19.7 + 5 log h) and SSRS2+CfA2 
130 (r < 130.0 h~ l Mpc, M B < -20.3 + 51og/i; |da Costa et al. 1994b| ), and one apparent- 
magnitude-limited sample CfA 101m (r < 101.05 h~ l Mpc, m^(o) < 15.5; Park et al. 1994 ). 



Now, while the Las Campanas galaxies have been selected from photometry in a red wave- 
length band, the SSRS2 and CfA2 galaxies were selected in the blue. As a rough guide to 
compare the red and blue magnitudes, we note that M^ CRS = —20.3 ( [Lin et al. 1996|) , while 
M£ fA2 = -18.8 ( Marzke et al. 1994| ) and M* SSRS2 = -19.5 ( |da Costa et al. 1994a|) . Nonethe- 
less, the LCRS P{k) is in remarkably good agreement with the P{k) of the magnitude-limited 
CfAlOlm and volume-limited SSRS2+CfA2 101 samples. There are some differences, though 
they are not at high significance given the errors. Specifically, the LCRS power spectrum 
is lower than SSRS2+CfA2 101 on intermediate scales A ~ 20 — 40 h^ 1 Mpc, and it is 
also low relative to these other two samples on large scales A ^ 150 h~ x Mpc. In contrast, 
the power spectrum of the volume-limited SSRS2+CfA2 130 sample is higher than that of 
the LCRS on all scales examined. Most simply, these P{k) differences may just be due to 
cosmic variance that arises from the different volumes being sampled in the three surveys. 
Alternatively, luminosity bias might help reconcile some of the differences we see, in the 
sense that intrinsically brighter galaxies are more strongly clustered. We note that the two 
volume-limited SSRS2+CfA2 samples are restricted to galaxies brighter than M^^ A2 and 
Mg SRS2 , whereas the LCRS sample contains galaxies between approximately M LCRS + 2 and 
M LCRS — 2. The magnitude-limited CfAlOlm sample does contain galaxies both brighter and 
fainter than M^f A2 , and seems the best match to the LCRS. Park et al. (1994) find evidence 
for luminosity bias in the CfA survey, and we explore this possibility for the LCRS in § |j. As 
another alternative, the differences between the two surveys may reflect intrinsic clustering 
differences of galaxies selected in different wavelength bands, red for the LCRS vs. blue for 
the other two surveys, analogous to the effect observed for CfA galaxies vs. infrared-selected 
IRAS galaxies (see below). The galaxy sampling and surface brightness selection criteria of 
the LCRS are unlikely to produce a systematic bias in our power spectrum estimates over 



the scales probed, given the results of the test done in § |2.2j and the agreement of the 50- 
and 112-fiber results in Figure [|(b). 

We next compare the LCRS P(k) to that found from two redshift surveys of IRAS 
galaxies: the 1.2 Jy survey (60 /im flux < 1.2 Jy; [Fisher et al. 1993] ) and the deeper but 



more sparsely sampled QDOT sample (60 jum flux < 0.6 Jy; |Keldman et al. 1994| ). As seen 
in Figure BL the IRAS 1.2 Jy P(k) is lower than that of the LCRS over all scales probed, 
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particularly for A < 30 h" 1 Mpc. Fisher et al. (1993) have shown that the 1.2 Jy P(k) is very 
similar in shape to the CfA 101 P{k), but that the overall clustering amplitude for IRAS 
galaxies is lower. In redshift space cr 8 (1.2 Jy) = 0.8, <7 8 (CfA101) = 1.1, while the LCRS is 
in between, with ^(LCRS) = 1.0. The clustering strength seen in the 1.2 Jy survey may 
reflect the inherently weaker clustering amplitude of the IRAS galaxy population, which 
tend to be spirals that avoid rich clusters (e.g. |Strauss et al. 1992j ). We have selected from 
the LCRS a sample of emission galaxies, with [Oil] A3727 equivalent widths W\ > 5 A, 
which should also tend to be spirals, as well as a sample of galaxies with W\ < 5 A, which 
should include more ellipticals (see Table [I]). We compare the relative convolved power 
spectra of the whole, emission, and non-emission LCRS samples in Figure ^. We find that 
for scales A ~ 10 — 100 h^ 1 Mpc, the emission sample is about 30% more weakly clustered 
than the full sample, while the non-emission sample is about 25% more strongly clustered 
than the full sample, and these results are roughly independent of scale. The same trends 
persist up to A ~ 300 h~ l Mpc, though the errors become larger for A ^ 100 h~ l Mpc (here 
again the errors have been estimated from the variation among the four LCRS subsamples 
divided by hemisphere and right ascension). Though the difference is not as great as for 
IRAS vs. LCRS galaxies, our results are consistent with an intrinsically weaker clustering 
strength of spiral and emission-line galaxies on scales A ^ 100 h~ l Mpc. Galaxies with early 
morphological types have been observed to be more strongly clustered than late types in a 
fair number of previous smaller samples, exploring generally smaller clustering scales (e.g., 
Davis fc Geller 1976 , Mo et al. 1992| , |Lovcclay et al. 1995| ). Our current results are thus in 



qualitative agreement, though the lack of one-to-one correspondence between morphologies 
and emission strengths vitiates a more quantitative comparison. 

We note next from Figure |8] that the LCRS and QDOT results are in somewhat better 
agreement, except that the QDOT slope is steeper, in the sense of having more relative power 
for scales A ~ 100 — 200 hT 1 Mpc compared to scales A ~ 50 h^ 1 Mpc. Feldman et al. (1993) 
also point out a turnover in the QDOT P(k) for A ^ 150 h^ 1 Mpc. The LCRS P(k) shows an 
apparent peak at roughly similar scales, A ~ 100 h~ x Mpc, but the constraint on the slope 
of the LCRS power spectrum on larger scales is not strong, as it can accommodate both 
declining to flat power spectra for A ^ 200 h^ 1 Mpc (§ |3.1| ). Note finally that the QDOT 
and 1.2 Jy results differ somewhat for A > 100 h~ x Mpc. Feldman et al. (1994) attribute 
this to sampling fluctuations rather than to luminosity bias among IRAS galaxies (see also 



Mann et al. 1995 ). The sampling fluctuation argument is also indicated by the analysis of 
Tadros & Efstathiou (1995), who have calculated the power spectrum of the combined 1.2 
Jy and QDOT samples. 

Overall, for scales A ^ 100 h^ 1 Mpc, despite some apparent differences, the power 
spectra from the various surveys are broadly consistent with each other given the large errors 
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common to all the surveys. On smaller scales, the power spectra of the different surveys show 
similar shapes, but there are differences in amplitude which may reflect intrinsic variations 
in clustering strength for galaxies of different luminosities or types. We will explore the 
dependence of clustering on luminosity for LCRS galaxies next. Nonetheless, it is noteworthy 
that the independent, optically-selected surveys LCRS and SSRS2+CfA2, which also provide 
the largest galaxy samples, do in fact agree well over the range of scales probed. 



4. Luminosity Bias in the LCRS 

As we saw in § [T2| and Figure 0, the SSRS2+CfA2 130 sample (M B < -20.3), which 
contains intrinsically brighter galaxies than the SSRS2+CfA2 101 sample (M B < —19.7), 
has a power spectrum amplitude ~ 40% greater than that of SSRS2+CfA2 101 (Park et al.| 
1994| |da Costa et al. 1994b| ). Motivated by this evidence of luminosity bias, here expressed 



as the variation of power spectrum amplitude with galaxy luminosity in the SSRS2 and 
CfA2 surveys, we would like to test for evidence of the same in the LCRS. Recall that 
these two samples are volume-limited, with galaxies brighter than Mq^ A2 = —18.8 and 
Mg SRS2 = —19.5, whereas our earlier magnitude-limited LCRS sample contains galaxies 
both brighter and fainter than M^ CRS = —20.3. Volume-limited samples, defined by both a 
lower and an upper bound on absolute magnitude, are a more natural choice over apparent- 
magnitude-limited samples for the purpose of examining how the power spectrum changes 
as a function of absolute luminosity. Also, for volume-limited samples the selection function 
is simpler to deal with. However, the 1.3-mag band in the isophotal magnitude limits of 
the 50-fiber data restricts drawing volume-limited samples to inconveniently small distance 
and absolute magnitude ranges, so we will confine the analysis to the 112-fiber data, which 
have wider 2.7-mag isophotal limits. Within the same redshift limits 10000 km s _1 < cz < 
45000 km s _1 as applied to the earlier apparent-magnitude-limited full LCRS sample, we 
form seven volume-limited samples, with absolute magnitude limits —18.5 > M > —19.5, 
-19.0 > M > -20.0, . . ., -21.5 > M > -22.5. Table g summarizes details of these samples. 
In the power spectrum calculation, we set the weight function for galaxy i to tu(rj) = Wi 
(see § O). 



One complication that we note from Table ^ is that the different volume-limited samples 
do not necessarily overlap much in space with one another, and one might wonder whether 
any differences seen among the samples may be due to sampling fluctuations rather than 
to actual luminosity bias. In order to better isolate the luminosity dependence, we adopt 
a procedure similar to that applied by Hamilton (1988) in examining the variation of the 
correlation function amplitude with galaxy luminosity in the CfAl survey. We divide each 
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volume-limited sample in half by absolute magnitude (i.e., the —19.0 > M > —20.0 would 
be split into two samples with -19.0 > M > -19.5 and -19.5 > M > -20.0), calculate the 
ratio of the power spectra of the brighter and fainter subsets (which are located in the same 
volume of space), and finally multiply (or divide) these ratios in successive volume-limited 
samples in order to determine the relative clustering amplitudes as a function of luminosity. 
We arbitrarily set the clustering amplitude of the —20.0 > M > —20.5 sample, which con- 
tains Ml CRS = —20.3, to unity. We apply this procedure to the whole-survey volume-limited 
samples, and estimate errors from the variation among volume-limited samples from the four 
subsets of the whole survey divided by hemisphere and right ascension (recall Figure 0(a)). 
We also average our results over four roughly equal bins in log A, as shown in Figure jlffi 
A = 5 - 15, 20 - 50, 60 - 150, and 200 - 400 h' 1 Mpc. We note first from the figure that 
results for the bin at 200 — 400 h^ 1 Mpc are quite noisy and the errors are large enough 
to be consistent with no variation of clustering amplitude with luminosity. We will focus 
attention on the three bins in the range from 5 — 150 h^ 1 Mpc, where the errors are not 
so large. We also find that the variation of clustering amplitude with absolute magnitude 
is fairly scale-independent over this range of A. We see that for —18.5 > M > —21.0 the 
clustering amplitude changes little with absolute magnitude. For galaxies brighter than 
M = —21.0, the clustering amplitude appears to rise with luminosity; in particular, for 
—21.5 > M > —22.0, the brightest bin with still reasonably small errors, the clustering 
amplitude is 1.6 ± 0.5 that at M*. The departures from no luminosity bias are relatively 
mild, except in the brightest bin, where unfortunately the errors are also large. Nevertheless, 
there does appear to be about a 50% stronger clustering for galaxies brighter than about 
M* — 1 relative to those fainter, in general agreement with the amount of luminosity bias 



seen in the CfA2 and SSRS2 surveys. The trends seen in Figure [Hj are similar to analogous 
results for scales <^ 10 h^ 1 Mpc from the much smaller CfAl survey ( [Hamilton 1988] ). Our 
results are also similar to recent findings from the large-volume, but l-in-20 sparse-sampled 



Stromlo-APM survey ([Loveday et al. 1996| , [Loveday et al. 1995]) , except that we do not see a 



much weaker clustering amplitude for sub-M* galaxies that Loveday et al. find. The differ- 
ence might be due the fact that we do not use galaxies fainter than about M* + 2 that were 
included in the analysis of the Stromlo-APM data. In addition, recent correlation function 
analysis of the SSRS2 sample QBenoist et al. 1996] ) indicates a luminosity bias trend similar 
to that seen in the LCRS. We note however that the LCRS provides advantages with its 
combination of sample size, sampling density, and depth over other redshift surveys, and 
likely provides the best current sample to examine the luminosity bias effect. 

The fact that the luminosity bias we see for super-M* galaxies persists to fairly large 
scales argues for a primordial rather than a local environmental effect as its root cause. 
Specifically, our results are qualitatively consistent with the biased galaxy formation scenario, 
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in which brighter galaxies form out of higher peaks in the underlying mass distribution and 
are therefore more strongly clustered, because such higher peaks are naturally more clustered 
in models like CDM with an initially Gaussian field of mass density fluctuations ( [Bardeeri 



et al. 19"8E ; White et al. 1987|) . Also, the large-scale clustering bias seen earlier for our 



emission and non-emission samples is also tied to some extent into the luminosity bias issue, 
as we have found that our non-emission sample is dominated by super-M* galaxies and our 
emission sample by sub-M* objects ([Lin et al. 1996|) . If the emission (~ morphology) and 



luminosity bias effects are independent, as evidenced in the Stromlo-APM survey QLoveday 



|et al. 19"95|) , then some primordial influence would also be suggested in relation to the 
galaxy star formation rate, for which [Oil] A3727 is a rough indicator ([Kennicutt 1992|) . In 
the present paper we have simply laid out our initial observations of these type-dependent 
clustering differences; we plan a more detailed examination of the the effects we have found 
and their quantitative implications for biased galaxy formation models. In addition, note 
that our present results have been limited to redshift-space, and it will be useful, though 
expectedly more difficult and noisy, to factor out the redshift distortion effects by examining 
the clustering in real space, which we also plan to attempt. 

We can however draw some brief conclusions in the context of the recent semi-analytic 
galaxy formation models of Kauffmann et al. (1995). In particular, their work indicates that 
the bias factor b should depend only weakly on scale (also |Wcinbcrg 1995 ), consistent with our 



observations. Also, their models predict that the amount of luminosity bias should decrease 
with increasing erg of the underlying mass power spectrum. Specifically, their Figure 6 
indicates that on scales ~ 20 h^ 1 Mpc, galaxies with M ^ M* + 1 are more biased relative 
to those with M ^ M* + 3, by factors of about 1.5 for admass = 0.4 and ^ 1.2 for a% )mass = 1. 



For comparison, we can convert the relative power values of Figure [10] into relative bias values 
(recall power oc b 2 ), and then make the appropriate sum over the LCRS luminosity function 
in order to determine the average bias over some absolute magnitude range. We find for 
scales 5 - 150 hr 1 Mpc that the bias of M ^ M* galaxies is < 1.1 times that of M ^ M* + 2 
galaxies (the stronger clustering of bright galaxies is more than compensated by their reduced 
numbers in the luminosity function). Our results are on the weak luminosity bias end of the 
range of Kauffmann et al. models, and thus imply that <Js,mass should be high, rs 1. As this 
is near the observed galaxy as , it also indicates that LCRS galaxies may be unbiased tracers 
of the mass. We return to this point below in § |5.2j. 
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5. Comparisons with Cosmological Models 

We now compare the LCRS results against the expectations of cosmological models. 
We begin with a brief comparison using two particular N-body simulations, which allow us 
to examine the model power spectra over both small and large scales. Then we proceed to 
widen the explorable parameter space by concentrating on fitting the LCRS power spectrum 
on large scales in the linear regime, and compare against the predictions of several classes 
of CDM-motivated models. In combination with COBE and other large scale structure con- 
straints, the large scale LCRS power spectrum will help us delineate the allowed parameter 
space in these cosmological models. 



5.1. Comparison to N-Body Simulations 

We compare the LCRS power spectrum to that of two N-body models, both kindly 



provided by Changbom Park: (1) the ODM model of § ^2, with f2 = 0.4, h = 0.5, and 
bias factor 6 = 1; and (2) the CDM1 model, which has Q = 1, h = 0.5, and b = 1.5. The 
CDM1 model was computed with a particle-mesh code on a 324 3 mesh, physical comoving 
volume (388.8 h' 1 Mpc) 3 , and contains 162 3 CDM particles and 1, 201, 320 biased "galaxy" 
particles, chosen by the biasing scheme of Park (1991). Both models are normalized so that 
the galaxy er 8 = 1. 

The convolved redshift-space power spectra of the LCRS and the N-body models are 
plotted in Figure O. The three power spectra agree well with each other for A < 20 h^ 1 Mpc. 



On intermediate scales A ~ 30 to 50 h~ x Mpc, the CDM1 model matches the LCRS results 
somewhat better than the ODM model does. For wavelengths A > 60 h~ l Mpc, the LCRS 
power falls in between the CDM1 and ODM curves, with the LCRS results closer to that 
of CDM1 on the largest scales A > 100 h^ 1 Mpc. Though both N-body models match the 
LCRS on small scales, neither model provides quite the right amount of large-scale power: 
not enough power in the case of CDM1 and too much in the case of ODM. 

To be more quantitative in our comparison, we can use a rank-sum test to see how 
probable it is to draw the LCRS power spectrum from the population of N-body mock- 
survey power spectra. We choose the rank-sum test as it is simpler to apply than a \ 2 test; 
we need not make any assumptions about either the distribution of power at each k or about 
the correlations between the power at different k. At each wavenumber k, we first assign 
each of the 30 ODM mock-survey plus 1 LCRS samples a rank in order of increasing 
power (the index i denotes the sample). Then, for each of the 31 samples we combine the 
ranks at different k by forming the sum Si = J2k Ri,k- Finally we assign each sample an 
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overall rank Ri )Sum in order of ascending S^. Doing this results in the LCRS receiving a rank 
RLCB,s,sum = 2, giving us a two-tailed confidence interval of 100 — 2 x (2/31 x 100) = 87% for 
rejecting the null hypothesis that the LCRS power spectrum is drawn from the population 
of ODM power spectra. This is not very high significance, and it will be further weakened by 
the fact that since we only have one ODM simulation box, we do not sample the full range 
of variation as we would given an entire ensemble of ODM simulation boxes. Nevertheless, 
the rank-sum test does give us a simple quantitative measure of how well the ODM model 
can match the LCRS data. Likewise, the rank-sum test applied to 30 CDM1 mock-survey 
plus 1 LCRS samples gives RLCRS,sum = 10, where we now assign R i k in order of decreasing 
power because the LCRS power is greater than that of most of the CDM1 mock surveys. The 
two-tailed confidence interval for rejecting CDM1 is then just 35%. Repeating the rank-sum 
tests, but now just focusing on large scales k < 0.1 h Mpc -1 , we find RLCRS,sum = 1 for both 
ODM and CDM1, so now the rejection probabilities are 94% in both cases. The match of 
the CDM1 model to the LCRS is thus appreciably worse when considering only large scales, 
while the match of the ODM model is little changed. Neither model provides a completely 
satisfactory match to the LCRS power spectrum over all the scales which we probe. The 
LCRS power spectrum lies between that of the two models, and as we will see in the next 
section, the LCRS power spectrum on large scales is better matched by a linear CDM model 
with QqH « 0.3. 



5.2. Fitting the Linear Power Spectrum and Comparison to COBE 

For any particular structure formation model, the best comparison procedure with ob- 



servations is to draw out mock redshift surveys as we did in § ^]2|. However, it will be useful 
to fit our observations to analytic models in order to explore more of cosmological parameter 
space than permitted by Park's two N-body models. To do a proper fit to a nonlinear, 
redshift- space power spectrum given a linear, real-space model requires us to account for the 
following: (1) the nonlinear evolution of the power spectrum; (2) the bias relating the galaxy 
fluctuations to the underlying mass fluctuations; and (3) the distortions caused by mapping 
from real to redshift space. Though these three effects may be modelled analytically (Pea- 



cock &c Dodds 1994 ; |Jain et al. 1995| ), we will take a simpler route and concentrate on fitting 



just those observations on large scales in, at least approximately, the linear regime. We thus 
sacrifice the information we have on small scales for the simplicity of not having to model 
the nonlinear effects (which perhaps are best dealt with via N-body simulations). 



To specify our assumptions, note the following relations connecting the galaxy power 
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spectrum and the underlying mass power spectrum: 

P galaxy, real— space^}^) fb\k)Prnass,real—space\k) (^5) 

Pgalaxy , redshift— spacei}^) fz\k)Pgalaxy,real—space\k) i (^6) 

where fb(k) is some function describing the bias of the galaxy power spectrum relative to 
that of the mass, and f z {k) describes the effect of redshift distortions caused by galaxy 
peculiar velocities. In general, the particular functional forms of fb and f z will depend on 
the parameters of the cosmological model. We make the simple assumption that over the 
scales we consider, fb{k) = b 2 = constant, where b is the bias factor. For f z (k), there are 
two main regimes. As derived by Kaiser (1987), on large scales in the linear regime, the 
redshift-space power spectrum is amplified by the infall of galaxies from low- to high-density 
regions, and 

2 n 4/7 1 n 8/7 

(More commonly in the literature the fractions 4/7 and 8/7 in the above equation are replaced 
by 0.6 and 1.2 respectively; but see Lightman fc Schechter 1990| .) On the other hand, on 



small scales in the nonlinear regime, the redshift-space power spectrum is suppressed by the 
smearing of the clustering pattern from small-scale galaxy peculiar velocities. For example, 
if we assume that galaxies have a small-scale Gaussian peculiar velocity field with a one- 
dimensional dispersion a v , then 



J y 1 2 ka v /H V ; 



( Peacock 1992 ). However, here we will focus on the linear regime and so assume that 



galaxy , redshift— space 



(k) 



4/7 , n 8/T 
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mass, real— space, linear 



(k) (29) 



holds as a useful approximation, that is, the redshift-space galaxy power spectrum is just a 
constant times the real-space mass power spectrum. 

Since the density fluctuations in the galaxy distribution should be less than unity in the 
linear regime, and since we observe that the galaxy fluctuations in 8 h~ x Mpc radius spheres 
og = 1, we need to restrict our fits at least to k < 2rr/(several x 10 h^ 1 Mpc). We will 
use k < 0.16 h Mpc" 1 , corresponding to A > 40 h^ 1 Mpc. Though this choice is somewhat 
arbitrary, it is motivated by our desire to extend the fit range beyond the flat part of the 
convolved power spectrum (to improve the handle on the power spectrum shape), but not so 
far that nonlinear effects become important. That ( p9"D approximately holds over our k range 
has justification. Theoretical work (e.g., [Kaufl'mann et al. 1995| ; |Weinberg 1995| ) indicates 
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that constant linear biasing does hold on large scales. There is also empirical evidence, from 
the comparison of real- and redshift-space galaxy power spectra from the survey compilation 
of Peacock & Dodds (1994), and in particular, from the comparison of various galaxy survey 
power spectra (including LCRS) against the mass power spectrum derived from peculiar 
velocity catalogs over the range 0.05 < k < 0.2 h Mpc -1 flKolatt fc Dekel 1995|) . We 



will thus fit various model linear power spectra (to be described below) to the full sample 
magnitude-limited LCRS convolved power spectrum for 0.016 < k < 0.16 h Mpc -1 . One 
complication is that because of correlations among the data points on large scales, a rigorous 
fit requires us to calculate and account for the covariance matrix of the errors (e.g., |Feldmari 



et al. 1994]) . However, for simplicity we take a more empirical approach and estimate the 
errors for the fit parameters from the variation among our usual four LCRS subsamples and 
from the N-body mock surveys. 

We focus on the case of inflationary, dark-matter dominated cosmological models, in 
which case we may write 

Pmass,real—space,linear(k) B mass k T (fc) , (30) 

where B mass is a normalization factor, k n is the inflation-produced primordial power spectrum 
of spectral index n, and T(k) is a transfer function that modifies the primordial spectrum. 
Note that the shape of the power spectrum is parameterized by the index n and by the 
form of T(k), which in turn depends on the specific matter content of the universe. Now, 
because we are fitting over a relatively small range in k, and because the convolved power 
spectrum tends to flatten out over the same k range, regardless of the true power spectrum 
shape, the LCRS results will not constrain the true shape very strongly. On the other 
hand, given a specific power spectrum shape, the LCRS results do give fairly stringent 
constraints on the power spectrum normalization. Using our earlier mock surveys and the 
usual procedure of dividing into four subsamples, we estimate a la fractional error of about 
10% in the normalization, obtained by fitting to the full-sample LCRS power spectrum 
over our desired k range. However, there is still freedom in the choice of bias parameter b 
(equation |27|) so that we cannot fix an absolute normalization from our galaxy data. The 



absolute mass normalization B mass can instead be determined independently using results 
from the COBE DMR experiment. We can then turn the problem around and use our derived 
galaxy normalization Blcrs to determine the value of b needed to reconcile the LCRS and 
COBE normalizations via the relation 



2 £lr/ 1 Qir/ 

1 + T7-T- + 



3 6 5 b 2 



Bmass ~ B LC RS ■ (31) 



We note that theoretical models suggest a bias b near unity or larger (e.g., [Kauflmann et al. 
1995|) . Empirically, we see within the LCRS that emission/non-emission samples have bias 
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values within ±15% that of the full sample (power oc 6 2 ), and that super-M* galaxies have 
bias ~ 1.3 times that of fainter galaxies. Likewise, Peacock & Dodds (1994) find a relative 
bias between optical and IRAS galaxies of b optica i : b IR As = 1-3, and Kolatt & Dekel (1995) 
find a range 0.77 — 1.21 in relative bias parameters for various surveys, including LCRS. In 
particular, the latter authors find b LCR s = fio 6 /(0-99±0.13) for k = 0.05-0.2 h Mpc^ 1 . For 
Qq = 1, this corroborates the earlier indication from our luminosity bias results that LCRS 
galaxies may have b ~ 1. In any event, it appears that 0.7 — 1.3 may represent a reasonable 
range in (absolute) bias values for galaxies, and the assumption that LCRS galaxies are 
roughly unbiased should serve as a useful guide to help constrain the models. 

To narrow down the cosmological parameter space to be examined, we will consider 
four simple variants of the basic flat, inflationary CDM model with Q = 1 and n = 1. 
These are: (1) ACDM, flat CDM models with non-zero cosmological constant Qa = 1 — ^o; 
(2) OCDM, open CDM models with fl < 1; (3) CHDM, flat models with both CDM and 
hot dark matter in the form of massive neutrinos; and (4) TCDM, flat CDM models with 
a tilted primordial spectrum of index n / 1. These models were motivated by efforts to 
address various deficiencies in matching the "standard" Qq = 1, h = 0.5 CDM model to 
observational constraints (e.g., see [Liddle et al. 199"6a| ,b,c and references therein). For all 



these models we need the basic CDM transfer function ( |Bardeen et al. 1986| ) 



TcDM(q) = ln(1 2 + 34 g 34g) I 1 + 3 - 89q + {16Aq)2 + {5Mq)3 + ( 6 - 71 <?) 4 ]~ 1/4 > ( 32 ) 

where q = k/hT and V is the so-called shape parameter. For no baryonic matter V is simply 
Q h, but more generally (e.g. [Liddle et al. 1996c[ ) 



r = fl hexp(— Q B — Q B /Q ) , (33) 

where Qb is the fraction of critical density in baryons. We will adopt the big-bang nucle- 
osynthesis value Qb = 0.015//i 2 (|Copi et al. 1995Q . As a parenthetical remark, we note that 
in previous LCRS papers (e.g., [Landy et al. 1996| ; |Lin 1995] ) we used a slightly different CDM 



transfer function ( |Bond fc Efstathiou 1984) ; [Ffstathiou et al. 1992| ) which has VLb = 0.03 and 



T = Qoh rather than the current definition (^). The difference in T is only 6% for Qq = 1. 
For our four cases we have 

T(k) = / T c DM ^ ACDM > 0CDM > TCDM (M) 

1 } \ T CDM (q)D(q,n u ) CHDM ' 1 } 

where D(q, Q u ) is a somewhat complicated function describing the effects of massive neutrinos 
(it is given in detail in |Pogosyan fc Starobinsky 1995| and [Liddle et al. 1996b]) . For our power 



spectrum definition, the necessary COBE normalizations may be expressed via 

B mass = (27r 2 c/H ) 3+n 5 2 H , (35) 
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where the four-year COBE DMR data ( Bennett et al. 1996 ) give flLiddlc ct al. 1996T5| ,c; Viana 
fc Liddle 1996| ; Bunn & White, in preparation) 



10 5 5 



11 



0.785-0.05 lnf2 



1.94 n 
x exp[— 0.95(n 



1) - 0.169(n - l) 2 ] ACDM, TCDM, CHDM(f2 = 1) 



1 95 f2-°- 35 -°' 191nn o 



OCDM 



(36) 

where the la fractional uncertainty in Sjj is 15% (quoted from |Liddle et al. 1996b| ). Also, to 
provide an additional handle on our parameter space, we consider the constraints provided 
by the observed abundance of galaxy clusters. The cluster results may be expressed in terms 
of the (linear theory) value of ag, mass as follows ( |Viana fc Liddle 1996| |Liddle et al. 1996"5| ,c): 

-(0.59-0.16^0+0. 06Og) 

(37) 



08,; 



0.60 Q Q y " u— ; ACDM, TCDM 

(0.36+0.31^0-0.28^) QQj^|^ 

CHDM 



0.60 Q 
[ 0.60 + 0.2^/3 



where the 95% confidence uncertainties are +32% and —24% (but enlarged by the factor 
^o.26iog 10 n for ACDM and by ^o.i7iog 10 n for OCDM; see |Viana k Liddle 1996j ). 



We consider each of our four classes of CDM models in turn. 

(1) ACDM. Here we find a shape parameter T = 0.3 ±0.1 (the error is one sdom among 
the four subsamples). We plot in Figure |12|(a) the resulting bias b as a function of Qq, where 
the error bands are approximate 2a values obtained by adding in quadrature the estimated 
2a fractional uncertainties in the LCRS (20%) and COBE (30%) normalizations B. Four 
sets of curves are shown, drawn for V = 0.3, 0.2, and for fixed Hubble constants h = 0.5 
and 0.8. (We note that a variety of measures of h, using Cepheids, supernovae, Tully-Fisher 
distances, surface brightness fluctuations, and other methods, indicate that h ^ 0.5, and 
likely of order 0.6 — 0.8; see e.g. |Jacoby et al. 1992] , [Freedman et al. 1994j , |Schmidt et al. 



1994 , |Riess et al. 1995 ). Recalling our earlier discussion, we also draw horizontal lines at 
b = 0.7 and 1.3 to roughly indicate the range of reasonable bias values. In Figure 0(b) 
we plot the corresponding COBE-normalized a^ mass values {2a error bands) and the 95% 
confidence limits from the cluster abundance constraints. We find that the preferred values 
T = 0.3 and 0.5 <^ h <^ 0.8 indicate 0.4 <^ Q ^ 0.7, but that these models are anti-biased, 
b ~ 0.5, and also violate the cluster abundance constraints. On the other hand, for V = 0.2, 
we do find models at Q w 0.4 — 0.5 and h w 0.5 which are unbiased and which also meet 
the cluster constraints. An Qq — 1, T — 0.3 model is unbiased but requires an unreasonably 
low h ~ 0.3 and is marginal with respect to cluster abundances. 

(2) OCDM. The V fit is the same as in the ACDM case. We plot the corresponding 
T = 0.3 and h = 0.5,0.8 curves in the b and a &mass vs. fl planes in Figure [13]. The most 



-24- 



dramatic difference in the open CDM case vs. the flat ACDM case arises from the much 
weaker variation with Q of the COBE normalization. The parameter space opens up a bit 
and for 0.4 Q <^ 0.6 there are models which satisfy T = 0.3 and 0.5 <^ h ^ 0.8, are 
unbiased, and match the cluster abundances. 



(3) CHDM. In Figure |Tj we plot our results, b and o" 8jmass vs. neutrino density VL V , for 
two cases h = 0.5, 0.8. The h = 0.8 results are strongly anti-biased, b pa 0.3, violate the 
cluster abundances, and therefore appear ruled out. On the other hand, the h = 0.5 models 
are mildly anti-biased, b ~ 0.8, and satisfy the cluster abundances for f2„ ^ 0.2. The LCRS 
fit further constrains Q u to 0.2 ±0.1 (the error is one sdom among the four subsamples) for 
the case h = 0.5. 

(4) TCDM. Here we plot b and as, mass vs. the spectral index n, again for the two cases 
h = 0.5 and 0.8, in Figure [15|. For h = 0.5 and n ~ 0.75 the models are unbiased and also 
satisfy cluster abundances. The LCRS fit gives n = 0.7 ± 0.2. For h = 0.8 we instead need 
n <^ 0.6 to meet the constraints, but this appears to violate the direct COBE DMR fit to 
the spectral index, n = 1.2 ± 0.3 (la; [Bennett et al. 19"96|) . 
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In summary, all the simple CDM variants we have considered can provide models that 
simultaneously satisfy the LCRS large scale power spectrum shape constraints, have Hubble 
constant 0.5 ^ h ^ 0.8, give an approximately unbiased LCRS galaxy distribution on large 
scales, and meet the cluster abundance constraints. For each of our four CDM classes, we 
plot in Figure O specific example linear power spectra that meet the above criteria: 



(38) 



The model linear power spectra are shown both unconvolved and convolved, and compared 
to the corresponding LCRS results. Note the similarity in the shapes of the convolved power 
spectra, which we are actually fitting. However, not surprisingly, the true unconvolved spec- 
tra do in fact differ more in shape, and they also show more power compared to the directly 
deconvolved LCRS power spectrum for A ^ 200 h^ 1 Mpc. Also shown are corresponding 
COBE la error boxes, shifted to redshift space using equation (^9]) and the appropriate bias 
value. 

The acceptable regions of our CDM model parameter space appear to be the following. 
Flat Qq — 1, h — 0.5 models are still viable, if we add massive neutrinos with Vt„ w 0.2 — 0.3 
or a spectral tilt n pa 0.7 — 0.8. Alternatively, flat f2o + = 1 CDM models with V pa 0.2 
and Qq pa 0.4 — 0.5 will also do, and similarly for open CDM models with Qq ~ 0.4 — 0.6 
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and T pa 0.3. We also note that these models can provide <Ts >mass values near unity, which 
was suggested in the earlier comparison of our luminosity bias results with the predictions of 
Kauffmann et al. (1995). Similarly, these models can also meet the constraint cr 8mass f2Q 6 ~ 
0.7 — 0.8 from peculiar velocity data ( [Kolatt fc Dekel 1995| ). Additionally, they also satisfy 
the 95% lower bound on the age of the universe, 12 Gyr, as determined from the ages of the 
oldest globular clusters QChaboyer et al. 1996| ). 

We have by no means attempted to exhaustively compare our results to all the possi- 
bilities (e.g., we can incorporate tilt in the open models, add gravity waves, and so on), but 
have simply explored some simple variations on the CDM theme, and have seen that there is 
already room for a number of such viable models. The next generation of very large redshift 
surveys, such as Sloan ( |Gunn fc Weinberg 1995| ) and 2DF ( [Ellis 199~3"D should provide us 
with improvements on the large-scale galaxy power spectrum ( |Vogeley 1995| ). In conjunc- 
tion with improvements in degree-scale microwave anisotropy experments (e.g., see reviews 
by (Scott et al. 1995| ; [Bond 1995| ) , the convergence of the large scale structure and microwave 
background observational fronts will hopefully give us a more definitive answer. On the 
other hand, for smaller scales A < 100 h~ l Mpc, P(k) is already quite well constrained in 
the LCRS. In particular, any viable model must match these small scale results, as we have 
seen earlier in the case of our two N-body models. The complication here is that galaxies 
exhibit type-dependent clustering differences (as is the case for our emission samples and for 
our super-M* galaxies), which will make it more difficult to interpret which types of galaxies 
are actually being simulated by the models, in as much as the process of galaxy formation is 
still not well understood, and there is freedom of choice in how one's "galaxies" are defined 
and represented in the simulations. The connection between the N-body models and the 
statistics of a galaxy sample will remain a weak link which will require a more thorough 
understanding of the types of galaxies selected by a particular survey and their relation to 
the underlying density fluctuations. 



6. Conclusions 

We have computed the power spectrum for a magnitude-limited sample of 19305 Las 
Campanas survey galaxies, with average redshift z — 0.1 and absolute magnitudes —22.5 < 
M < -18.5, over scales A = 5 - 400 h' 1 Mpc. P(k) may be approximated by a power law 
P{k) oc k n , with slopes n = — 1.8 ± 0.1 on small scales A = 5 — 30 h~ l Mpc, and n = 1 ± 1 on 
the largest scales A ~ 200 — 400 h~ x Mpc. The change in slope of P(k) is real and could be 
the result of a possible maximum in P(k) observed at A ~ 100 h~ l Mpc, though our errors 
preclude a definitive detection of the turnover. The overall amplitude may be expressed as 
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a$ = 1.0 ±0.1 in redshift space. We have also compared our results to those from other 
redshift surveys. For A <^ 100 Yc 1 Mpc, the LCRS P(k) is similar in shape to that of the 
other samples, but with some differences in amplitude that may reflect intrinsic clustering 
differences for galaxies of different luminosities or types (below). For A ^ 100 h' 1 Mpc, the 
LCRS results are consistent with those of other surveys given the large errors among all the 
surveys on these scales. The overall LCRS results agree best with that from the combined 
SSRS2+CfA2 redshift surveys. 

We find evidence in the LCRS for the dependence of clustering strength on galaxy type. 
A sample of LCRS emission galaxies with [Oil] A3727 equivalent widths > 5 A are about 
30% less clustered than the full LCRS sample on scales A <^ 100 h^ 1 Mpc. This may reflect 
an intrinsically weaker clustering strength for spiral galaxies, in qualitative agreement with 
the weaker clustering of IRAS vs. optically-selected galaxies. We also compute the power 
spectra of volume-limited LCRS samples to test for evidence of luminosity bias. We find that 
LCRS galaxies in the range —18.5 > M — 5 log h > —21 appear to be clustered similarly, over 
scales A ~ 5 — 150 h^ 1 Mpc. However, more luminous galaxies, which lie about 1 magnitude 
brighter than M^ CRS = —20.3 + 5\ogh, exhibit about a 50% stronger clustering amplitude, 
in agreement with similar results seen in the CfA2 and SSRS2 surveys and the qualitative 
expectations of biased galaxy formation scenarios. 

Our full nonlinear power spectrum is intermediate between that of two CDM N-body 
models, an open Q h = 0.2 model and a standard QqH = 0.5 model, both normalized to a 8 = 
1 for galaxies. The models agree well with the LCRS results on small scales A ^ 50 h~ l Mpc, 
but neither model provides a completely satisfactory fit to the full LCRS power spetrum. 
To explore more of model parameter space, we then fit the large scale LCRS power, on 
approximately linear scales A ^ 40 h' 1 Mpc, to several classes of linear CDM models. We 
find that a number of viable models remain which can simultaneously satisfy constraints 
provided by the LCRS power spectrum, the Hubble constant range 0.5 ^ h ^ 0.8, cluster 
abundance results, and the reasonable assumption that LCRS galaxies are approximately 
unbiased on large scales, relative to the mass normalization provided by the 4-year COBE 
DMR data. These models include: (1) flat fi A + fi = 1 CDM models with fi « 0.4 - 0.5 
and shape parameter T « Qoh 0.2; (2) open CDM models with Qq rs 0.4 — 0.6 and T m 0.3; 
(3) flat Qo — 1 models with CDM plus massive neutrinos of density Vt v ~ 0.2 — 0.3; and (4) 
flat tt = 1 CDM models with a spectral tilt n « 0.7 - 0.8. 

The present situation for power spectrum measurements from redshift surveys appears to 
fall into two regimes. On scales A <^ 100 h~ x Mpc, the power spectrum is well determined, and 
various surveys have consistently detected differences in the clustering strengths of different 
galaxy populations, though the latter effect may complicate interpretation of model predic- 



-27- 



tions. On the largest scales A ^ 100 h~ x Mpc, however, the errors in P(k) are still large. We 
have not been able to conclusively rule out any of the popular classes of CDM models, though 
we have begun to constrain their parameter space. Nonetheless, the observational picture 
should improve when the next generation of very large redshift surveys come on line, and 
as the degree-scale microwave background observations progress. In the meanwhile, the Las 
Campanas results have already provided constraints on P(k) on both large and small scales, 
complementary to, and largely independent of, earlier results for smaller and shallower sam- 
ples. Furthermore, the type-dependent clustering differences observed in the survey should 
provide useful observational constraints on models that seek to reproduce accurately the 
clustering properties and formation histories for a realistic mix of galaxy types. Toward this 
end, interested readers are invited to obtain the Las Campanas Redshift Survey catalog, now 
publicly available at the Internet site "|http: / / manaslu.astro.utoronto.ca/ "lin/lcrs-htmlf . 



We thank Changbom Park for providing the N-body models, Michael Vogeley for useful 
discussions and for providing the SSRS2+CfA2 results, Hume Feldman for the QDOT results, 
and Tsafrir Kolatt for helpful comments on the manuscript. The Las Campanas Redshift 
Survey has been supported by NSF grants AST 87-17207, AST 89-21326, and AST 92-20460. 
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Note Added in Proof: Power spectrum results for the Stromlo-APM redshift survey 
( [ladros fc Efstathiou 1996]) appeared after submission of the present paper. Time does not 
permit a detailed comparison, but we note that the Stromlo-APM authors report that their 
results are consistent with those of the CfA2 101 sample ( [Park et al. 19941 ) and that they 
also find evidence for the stronger clustering of galaxies brighter than M*. 
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Table 1. Magnitude-Limited Samples 



Sample 


N qa l 


All 


19305 


North RA < 12 h 7 


5033 


North RA > 12 h 7 


4327 


South RA < h 7 


5125 


South RA > h 7 


4820 


All 50-fiber 


4070 


All 112-fiber 


15235 


3727 W x > 5 A 


7712 


3727 W\ < 5 A 


11593 



a All samples further restricted to 10000 km s -1 < cz < 45000 km s _1 
and -22.5 < M - 51og/i < -18.5 ("hybrid" Kron-Cousins R magni- 
tudes). 



Table 2. Volume-Limited Samples 



Absolute Magnitude Limits a 


Distance Limits b (/i 1 Mpc) 




-18.5 > M > -19.5 


98 < r < 160 


732 


-19.0 > M > -20.0 


98 < r < 198 


1867 


-19.5 > M > -20.5 


119 < r < 243 


3179 


-20.0 > M > -21.0 


147 < r < 298 


4591 


-20.5 > M > -21.5 


181 < r < 362 


4404 


-21.0 > M > -22.0 


224 < r < 405 


2534 


-21.5 > M > -22.5 


274 < r < 405 


582 



a "Hybrid" Kron-Cousins R magnitudes with h—1. 

b For nominal 112-fiber apparent magnitude limits 15.0 < m < 17.7. 
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Fig. 1. — The LCRS galaxy distribution in the Northern and Southern galactic caps. 
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Fig. 2. — The true unconvolved and convolved power spectra of the ODM N-body model, and 
the corresponding power spectra as measured from the ODM LCRS mock surveys. The error 
bars are one standard deviation of the mean (sdom) as measured from the mock surveys. 
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Fig. 3. — Integrand k' P(k')K(k, k') of the convolution integral (16) for the ODM model at 
5 values of k. Note we have multiplied the integrand by A;' In 10 A(log 10 /c') as appropriate 
for our log-linear plot, and we have also normalized by dividing by P(k). 
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Fig. 4. — (a) The directly measured convolved power spectrum for the full magnitude-limited 
LCRS sample and for each of the four subsets of the full sample divided by hemisphere and 
right ascension, (b) The convolved power spectrum for the 112- and 50-fiber subsets of the 
full sample. 
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Fig. 5. — The effects on the power spectrum of the full magnitude-limited LCRS sample 
arising from: (a) assigning fake velocities to unobserved objects lying within 55" of an 
observed galaxy; and (b) applying geometric weights to correct for the reduced spectroscopic 
success rate at field corners. See text for details. 
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Fig. 6. — The fit to the LCRS power spectrum using the model of equation (23). 
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Fig. 7. — The deconvolved LCRS power spectrum and fit compared to the power spectra of 
three samples drawn from the combined SSRS2 and CfA2 redshift surveys. 
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Fig. 8. — The deconvolved LCRS power spectrum and fit compared to the power spectra of 
the IRAS 1.2 Jy and IRAS QDOT redshift surveys. 
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Fig. 9. — The convolved power spectra of two LCRS subsamples, defined by [Oil] A3727 
equivalent width W\, relative to that of the full sample. The errors shown are standard 
deviations of the mean determined from the four LCRS subsets divided by hemisphere and 
right ascension. 
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Fig. 10. — The relative convolved power as a function of absolute magnitude, normalized to 
one at M — —20.25, shown for the full LCRS sample (heavy solid lines) and for the four 
LCRS subsamples divided by hemisphere and right ascension (light broken lines). The errors 
shown are one standard deviation of the mean from the four subsamples. 
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Fig. 11. — The observed convolved LCRS power spectrum compared to the convolved power 
spectra of the N-body models ODM and CDM1. See text for details. 
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Fig. 12. — (a) Bias b vs. f2 for ACDM models where Q + Q A — 1. 2a error bands are 
shown for fixed shape parameters T and fixed Hubble constants h. The lines b = 0.7 and 1.3 
roughly indicate the range of reasonable bias values, (b) COBE-normalized <J8 tmass values vs. 
Qq for the same models as in (a). The unlabeled thin dashed lines indicate 95% confidence 
limits from cluster abundance constraints. See text for more details. 




Fig. 13. — Similar to Figure [12| except that open CDM models with fi < 1 are being 
considered. 
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Fig. 14. — Similar to Figure [12] except that f2 = 1 CHDM models with massive neutrinos 
of density Q u are being considered. The vertical lines indicate the LCRS best fit range 

0.5. 



Q v = 0.2 ± 0.1 (la) for the case h 
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Fig. 15. — Similar to Figure |12| except that Q = 1 tilted CDM models with spectral index 
n are being considered. The vertical lines indicate the LCRS best fit range n = 0.7 ± 0.2 
(la) for the case h = 0.5. 
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Fig. 16. — The LCRS power spectrum and the fits to it using the linear CDM models 
described in the text, shown both in convolved and unconvolved form. Also shown are 
corresponding COBE la error boxes, transformed into redshift space using the indicated 
bias values. See text for details. 



